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Abstract 

We study the dynamics of a straight vortex line in a partially Bose-Einstein condensed atomic 
gas. Using a variational approach to the stochastic field equation that describes the dynamics 
of the condensate at nonzero temperature, we derive the stochastic equations of motion for the 
position of the vortex core. Using these results, we calculate the time it takes the vortex to spiral 
out of the condensate. Due to the fact that we include thermal fluctuations in our description, this 
lifetime of the vortex is finite even if its initial position is in the center of the condensate. 
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I. INTRODUCTION 



Contrary to classical fluids, superfluids support rotation only through quantized vortices. 
Since quantized vortices are therefore one of the hallmarks of superfluidity, the experimental 
and theoretical study of these topological excitations has attracted a great deal of attention 
in the field of Bose-Einstein condensed gases. Following their first experimental observation 
they have now been observed and experimentally studied by various groups @,IEIH]- 

Theoretically, the dynamics of a single vortex line in a Bose-Einstein condensate has been 
studied extensively in the zero-temperature limit jE}. In the absence of external rotation, 
the vortex is predicted to precess around the center of the condensate, which has indeed 
been observed experimentally by Anderson et al. p]. However, in this experiment it is also 
observed that the distance of the vortex core to the center of the condensate increases with 
time, i.e., the vortex spirals out of the condensate. This observation is a sign of the presence 
of dissipation and thus cannot be explained on the basis of a zero-temperature approach. 
To understand it, we have to include the effects of the noncondensed component of the gas. 
Although the effects of the noncondensed thermal cloud on the equilibrium properties of 
rotating Bose gases have been investigated 0, H, H, EIJ , the nonequilibrium dynamics of a 
single vortex at nonzero temperatures has attracted relatively little attention. 

It is the purpose of this paper to study the effects of the thermal cloud on the motion 
of a single vortex in a Bose-Einstein condensate. The star ting point of our study is the 
stochastic Gross-Pit aevskii equation derived by one of us |ll|, ll2|. This equation generalizes 
the usual Gross-Pitaevskii equation, that provides an accurate description of the dynam- 
ics of the condensate at zero temperature, to nonzero temperatures by the inclusion of a 
dissipative term that describes the growth or decay of the condensate. Moreover, thermal 
fluctuations are included by an additive noise term that is related to the growth or decay of 
the condensate by means of a fluctuation-dissipation theorem. Although this equation can 
be solved numerically we intend to capture as much of the physics as possible by using 
a variational approach to stochastic field equations that was developed by two of us "Lll . 



To make the paper more self-contained, we briefly discuss the zero-temperature dynamics 
of a vortex in a Bose-Einstein condensate in SeclHI In Sec. Illll we then derive the stochastic 
equations of motion that describe the motion of the vortex core at a nonzero temperature. In 
Sec. |TV|we use these equations to derive an equation for the average distance of the vortex to 
the center of the condensate and use the latter result to calculate the lifetime of the vortex. 
We compare our results with the available theoretical results 2> 3- Unfortunately, there 



does not exist a detailed experimental study of the lifetime of the vortex. We end in Sec. |V| 
with our conclusions. 



II. VORTEX DYNAMICS AT ZERO TEMPERATURE 



The dynamics of a Bose-Einstein condensate is, at zero temperature, well described by 
the Gross-Pitaevskii equation 
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for its macroscopic wavefunction 0(x, t). Here, m is the atomic mass, and T 2B = Anah 2 /m 
is the two-body T(ransition) matrix at zero energy, with a > the s-wave scattering length 
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of the atoms of interest. We take the external trapping potential of the form 



y ext (x) = -m 
v ; 2 
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with u) z ^> u, which implies that we have a pancake-shaped condensate. The reason for 
choosing this geometry is that it allows us to neglect the curvature of the vortex line, 
provided that the vortex is close to the center of the condensate. 

To study the dynamics of a vortex at zero temperature we use the variational ansatz 
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where N c is the number of condensate atoms. This ansatz has the same form as the exact 
wave function for a noninteracting Bose-Einstein condensate with a vortex along the z-axis. 
We treat the width of the condensate in the radial and axial directions, denoted by q and 
q z , respectively, as time-independent variational parameters. The pancake-shaped geometry 
implies that q z <C q. The coordinates of the vortex in the x-y plane, denoted by Ui(t), 
are taken to be time-dependent variational parameters. For the case that the vortex is 
close to the center of the condensate, we expect that the variational ansatz in Eq. © offers 
a quantitatively good description of the wave function of the condensate in the weakly- 
interacting limit that is determined by the condition 
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with N c the number of atoms in the condensate. 

To determine the equations for the variational parameters we note that the Gross- 
Pitaevskii equation follows from a variation of the action 
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After substitution of the ansatz in Eq. into this action we find the result 
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with the precession frequency of the vortex equal to 
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and the potential, which represents the total energy of the condensate per particle, given by 



h 2 h 2 o 9 1 9 aN c h 2 

V{q, q z , N c ) = — + — — + muj 2 q 2 + -muj z q 2 z + - _ „ . (8) 



mq 2 Amq 2 4 2\/2iimq 2 q 

As a result we observe that the equilibrium widths of the condensate in radial and axial 
direction, denoted from now on again by by q and q z , respectively, are determined by min- 
imizing the potential V(q, q z , N c ) for a given number of condensate atoms. By varying the 
action with respect to the position of the vortex we find that the motion of the vortex is 
determined by the equations 

u x (t) = -uj p (N c )u y (t) , 

u y {t) = u p {N c )u x {t) , (9) 

which imply that the vortex precesses around the center of the condensate with precession 
frequency u p (N c ), which is equal to u p (q,q z , N c ) evaluated at the equilibrium values of q 
and q z . 

Although the above variational analysis has already provided us with a simple description 
of the precession of the vortex around the center of the condensate, it fails to account for 
the experimental obervation that the vortex spirals out of the condensate in the absence of 
rotation |6]. Contrary, the vortex remains, according the equations of motion in Eq. Q, 
at a fixed distance from the condensate center. In particular, this implies that if the initial 
position of the vortex is the center of the condensate, it will remain there forever. This 
discrepancy between theory and experiment is resolved by including the effects of thermal 
fluctuations, as we will see in the next sections. 



III. VORTEX DYNAMICS AT NONZERO TEMPERATURE 

The dynamics of a Bose-Einstein condensate is at nonzero temperatures, i.e., in the 
presence of a thermal cloud at temperature T = l/(kn0) and with a chemical potential /i, 
determined by the stochastic Gross-Pitaevskii equation [n, 2, k| 3 



x|-^ + V ox, (x)-,. + r 2B Wx,i)| 2 U(x,«) + ^x,i) . (10) 

This Langevin field equation quite generally generalizes the Gross-Pitaevskii equation to 
nonzero temperatures, and includes both dissipation, i.e., decay or growth of the Bose- 
Einstein condensate, and (thermal) fluctuations. The dissipation and the noise arise physi- 
cally due to collisions between noncondensed atoms in which one of the atoms is scattered 
into the condensate, and the time-reversed process. The complex gaussian noise in the 
stochastic Gross-Pitaevskii equation is completely determined by the correlations 

( V *(x,t) V (x',t')) = ^S K (x,t)5(t-t')5(x-x') , (11) 
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where the strength of the noise is determined by the so-called Keldysh self-energy, given by 
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In this expression, is the Bose-distribution function of the thermal cloud, evaluated at 
an energy of a thermal particle, which is in the Hartree-Fock approximation given by 
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This approximation is valid in the regime where the temperature is large compared to the 
zero-point energy of the external trapping potential but smaller than the critical tempera- 
ture. Note also that we have taken the thermal cloud to be in equilibrium, although this 
can easily be generalized. 

In first approximation we neglect the inhomogeneity of the thermal cloud over the size of 
the condensate. In good approximation, the Keldysh self-energy then turns out to be given 

I'.v 

which we will use from now on. Note that since the Keldysh self-energy determines both 
the dissipation and the strength of the fluctuations, the stochastic Gross-Pitaevskii equation 
can be show to automatically fulfill the fluctuation-dissipation theorem which ensures that 
the condensate relaxes to the correct equilibrium distribution. 

The (Wigner) probability distribution for the condensate wave function, that results from 
the stochastic field equation in Eq. (|1U|). can be written as the functional integral 



with the effective action given by 
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To determine the dynamics of the vortex at nonzero temperatures, we have to slightly 
generalize the ansatz in Eq. (jSJ) to allow also for fluctuations in the number of condensate 
atoms and the global phase of the condensate wave function Hence we now use the 

ansatz 
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For the moment we consider the case T 2B = 0, but keep the Keldysh self-energy nonzero. 
This implies that we are dealing with a noninteracting Bose gas in contact with a heat bath. 
Subsitution of the above ansatz into the action leads to 
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where we have again taken for the time-independent variational parameters q and q z the 
values obtained by minimizing the potential in Eq. (JHJ). The damping rate 7 is given by 



2„ T K 
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and the condensate chemical potential reads 
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From the action for the variational parameters we are able to deduce the stochastic rate 
equation for the number of atoms and the stochastic equations of motion for the global 
phase and for the position of the vortex. First, the stochastic rate equation for the number 
of atoms is given by 



dNJf) 
dt 



-^£ K [fi c (t) - fi] N c (t) + 2^Ut) V {t) 



with the correlations of the Gaussian noise given by 
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Second, the stochastic equation of motion for 9o(t) is given by 14] 
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where the correlation of the gaussian noise are given by 
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(23) 
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Since we are mostly interested in the stochastic equations of motion for the position of 
the vortex, we neglect from now on the fluctuations in the number of condensate atoms and 
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the global phase of the condensate. The Langevin equations for the position of the vortex 
core are, from the action in Eq. ([18)1 . seen to be given by 



u x (t) = -ujp(0)u y (t) + ju x (t) + T) x (t) , 

u y (t) = uj p (0)u x (t) + JUy(t) + T]y(t) , (25) 

where the gaussian noise is completely determined by 

(ViWViV)) = if" - = 'Wit - • (26) 

The stochastic equations for the position of the vortex core in Eq. (|23|) have two extra terms 
with respect to the zero-temperature result in Eq. (jHJ). First, the stochastic equations have 
a damping term proportional to the damping rate 7. As we will see, this term leads to 
exponentional increase in the distance between the vortex position and the center of the 
condensate. Secondly, the stochastic equations have an additive noise term which represents 
the effect of thermal fluctuations on the position of the vortex core. If these noise terms 
would not have been present, the lifetime of the vortex would again be infinite if the initial 
position of the vortex is the center of the condensate. Physically, the vortex is "kicked" out 
of the center of the condensate due to thermal fluctuations which are represented by the 
gaussian noise terms in Eq. (J25|) . 

Although we have in first instance derived the stochastic equation for the position of the 
vortex core by assuming that we are dealing with a Bose-Einstein condensate of noninteract- 
ing atoms, we can now consider the interacting case by replacing the precession frequency in 
cUp(O) by lu p (N c ). The Fokker-Planck equation, that determines the probability distribution 
for the position of the vortex, is then given by 
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In the next section, we use this equation to determine the average distance of the vortex to 
the center of the condensate. This result is then used to obtain the lifetime of the vortex. 



IV. LIFETIME OF THE VORTEX 

The average distance of the vortex to the center of the condensate is given by 

r (f) = (\/ u l + u l) (*) = J du x du y ^u 2 x + u 2 y P[u x , u y , t] . (28) 

With the Fokker-Planck equation for the position of the vortex core in Eq. ()27jl we find, 
with the use of a partial integration, that r(t) obeys the equation of motion 
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The general solution of this equation is given by 



/[2 7 r2(0) + q]e^-q 
r (t) = 7^= • (30) 



Note that for small times t and distance r(t) we have that r(t) oc v^- This may be understood 
from the fact that, if we neglect the contribution of the precession frequency u> p (N c ) and 
the damping rate 7 to the Fokker-Planck equation in Eq. ()27|). the vortex simply undergoes 
Brownian motion. For large times t and distance r(t) we have that r(t) oc e 7 *, as is easily 
seen from Eqs. fj2Uj) and (JHUj) . 

For a given initial position r(0) = r min , the vortex lifetime is defined as the time it takes 
the vortex to reach the edge of the condensate at r max . Of course, for our gaussian ansatz the 
edge of the condensate is not defined uniquely, but we could, for instance, take r max = q. In 
the Thomas-Fermi limit we have r max = R, with R the Thomas-Fermi radius. The lifetime 
of the vortex is in first instance given by 

'r 2 max + a/(2 1 y 



-im 

27 
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Let us first discuss the case where we neglect the contribution of the noise due to thermal 
fluctuations, i.e., we take a = 0. The lifetime of the vortex is then given by 

T^it'-YTAr^C—) • 02) 



6 \aj (k^,T) 

This result is very similar to the result obtained bv Fedichev and Shlyapnikov [15] who in- 
deed neglect thermal fluctuations. Schmidt et al. (l6|, who include fluctuations, numerically 
find a similar result. In particular, the dependence on the ratio r max /r min is identical, and 
we observe that, in the absence of thermal fluctuations, the lifetime of the vortex diverges 
as r min — > 0. Interestingly, although we consider the weakly-interacting limit, as opposed to 
Fedichev and Schlyapnikov who consider the Thomas- Fermi limit, the temperature depen- 
dence of the prefactor of the lifetime in Eq. ()32|) is the same in both results. This is even 
more surprising because of the fact that the dissipation of the vortex considered by Fedichev 
and Shlyapnikov is physically due to scattering of quasiparticles of the vortex core. This 
mechanism of decay is physically quite distinct from to the collisional effects we consider 
here. 

Including thermal fluctuations, the vortex lifetime is finite even if its initial position is 
the center of the condensate. In particular, for that case we have that 



7r / q\ 2 h 
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We observe that the leading low-temperature dependence of the lifetime of the vortex is 
T _1 In [N c hu /(kBT)], whereas in the case where we neglected the thermal fluctuations we 
had that the lifetime was proportional to T~ x . 



V. CONCLUSIONS 

We have studied, by means of a variational analysis, the motion of a vortex in a Bose- 
Einstein condensate. The main results were the stochastic equations of motion for the 
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position of the vortex and their corresponding Fokker-Planck equation, which enabled us to 
analytically calculate the lifetime of the vortex. We have compared our results for the latter 
with the available theoretical results 0,0]. It is surprising that, although the collisional 
decay mechanism that is considered here is physically quite different from the decay mecha- 
nism discussed by Fedichev and Shlyapnikov [15( , we nevertheless find the same temperature 
dependence if we neglect the thermal fluctuations. Including thermal fluctuations leads to 
a different temperature dependence. Note also that, because we consider condensate dis- 
sipation due to collisions, the lifetime of the vortex is inversely proportional to the square 
of the scattering length. The result for the lifetime found by Fedichev and Schlyapnikov is 
proportional to a -4 / 5 . It would therefore be very interesting to experimenally measure the 
dependence on the temperature and the scattering length of the vortex lifetime. 

Our present analysis was limited to the case where the overall density profile is close 
to a gaussian, i.e., the so-called weakly-interacting limit. In principle, our treatment can 
be generalized to the case where the condensate is described by a Thomas-Fermi density 
profile. However, in this limit technical problems occur due to the fact that we need to 
incorporate the density profile of the core of the vortex. This strongly limits the feasibility 
of analytical results. We expect, however, that the main difference between the weakly- 
interacting limit and the Thomas-Fermi limit is the difference in the precession frequency 
of the vortex, which can easily be found from a zero-temperature approach 0. Therefore, 
we believe that the lifetime results presented here provide also in the Thomas-Fermi limit a 
more than qualitative understanding of the dynamics of a vortex at nonzero temperatures. 
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